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Abstract In this study, non-radioactive Digoxigenin labeled ribosomal DNA (rDNA) probes 
were used for Southern blotting analysis to study the molecular phylogeny of the giant panda and related 
species. Restriction maps in the regions of rDNA spacers were compared between giant panda ( Ail- 
uropoda melanoleuca }, lesser panda (Ailurus fulgens), Asiatic black bear ( Selenarctos thibetanus ) , 
sun bear ( Helarctos malayanus}, raccoon (Procyon lotor ) and lynx ( Felis Lynx}. Phylogenetic trees 
for these species were constructed using maximum likelihood and parsimony method. The results show 
that in respect to (DNA RFLPs, the giant panda is more closely related to bear than to lesser panda; 
while the lesser panda is slightly related to the raccoon. 
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1 Introduction 


The giant panda {Ailuropoda melanoleuca ) is an endemic species of China and might be 
the most popular wild animal in the world. As matched with its popularity, numerous studies 
have been done on the morphology, ecology and behavior of the giant panda. However, al- 
though the giant panda is regarded as a “living fossil” in zoology and paleontology, its evolu- 
tionary relationships with other related species in Arctoidea, such as with the bear or the lesser 
panda, have been a biological puzzle (Van Valen, 1986). Most classification studies vary in 
placing the giant panda as a bear, a relative of the lesser panda, or a distinct family (Van 
Valen, 1986; Zhang et al., 1993). 

In the past decades, traits related to anatomy, behavior, karyology, immunological dis- 
tances, isozyme electrophoresis and DNA hybridization, were used to clarify the classification 
and evolution of the giant panda (Sarich, 1973; O’Brien et al., 1985; Van Valen, 1986). 
In 1985, O’Brien et al. analyzed the molecular genetic data obtained from four different ge- 
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netic methods and seemed to have solved the dilemma by concluding that the giant panda is 
closer to the bear. Subsequently, research with hemoglobin sequences (Tangle et al., 1986) 
and mitochondrial DNA (mtDNA) RFLPs (Zhang et ai., 1991), placed the giant panda as a 
close relative of the lesser panda. More recent analyses on mtDNA sequences (Zhang et al., 
1993) as well as maximum likelihood analysis of e-hemoglobin and B-hemoglobin sequences 
(Hashimoto et al ., 1993) have again showed that the giant panda is closer to bear. More da- 
ta, especially from nuclear DNA, would aid in consolidating the “giant panda-bear” relation- 
ship. To our knowledge, the nuclear DNA RFLPs of the giant panda and related species have 
not been studied yet. 

Ribosomal RNA (rRNA) genes are multigene families, and generally there are several 
hundred copies in one mammalian genome (Amheim et al., 1983). The copies are usually 
distributed in tandem repeats at several different sites on chromosomes. Each rDNA repeating 
unit is composed of three rRNA genes, 28S, 5.8S and 18S RNA, which are separated from 
each other by spacers. The spacer regions, especially the non-transcriptional spacer between 
28S and 18S RNA genes, are known to evolve rapidly, and considerable RFLPs are apparent 
between populations and species (Arnheim et al., 1983; Wilson et al ., 1984; Hillis et al., 
1986, 1988; Suzuki ez až., 1990, 1994). Because of the strong forces that encourage homog- 
enization within a genome as well as within a population, sequence homogeneity is maintained 
within a population (Dover, 1982). Therefore, small samples could usually be a good repre- 
sentative in phylogenetic analysis of a species or a population (Hillis e ai., 1986). These fea- 
tures described above make the analysis of (DNA RFLPs a good approach for phylogenetic stud- 
jes at the species, genus and even family level, and a lot of achievements have been reached 
(Arnheim et al., 1980; Wilson et al., 1984; Hillis et af., 1986, 1988; Suzuki et al., 
1990, 1994; Allard et al., 1991; Wall ez af., 1992; Lan er al., 1993). 

In this study we compared the restriction maps for 13 restriction enzymes on the rDNA 
spacer regions of the giant panda and related species using an non-radioactive DNA labeling and 
detecting method (Boehringer Mannheim). Phylogenetic reconstruction using maximum likeli- 
hood and maximum parsimony were conducted with the resultant matrix of sites. The results 
suggest that phylogenetic relationships among giant panda, lesser panda, bear and raccoon are 
similar to those based on mtDNA and hemoglobin sequence comparisons (Zhang et al., 1993; 
Hashimoto et al., 1993). 


2 Materials and Methods 


2.1 Samples 

Blood samples of the giant panda Aifuropoda melanoleuca and the raccoon Procyon lotor 
were provided by the Chengdu Zoo, China. Samples of the other species, lesser panda Ailurus 
fuigens, Asiatic black bear Selenarctos thibetanus, sun bear Helarctos malayanus and lynx 
Felis Lynx, were collections of the Kunming Cell Bank, the Chinese Academy of Sciences. We 
used two individuals of giant panda and lesser panda in this study, and one for the remaining 
species. 
2.2 Southern hybridization analysis 

DNA was isolated from bloods or cultured fibroblasts according to the method described by 
Sambrook er al. (1989). Genomic DNA (2-5 mg) was digested with 5 U restriction en- 
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donuclease at 37°C overnight. Completely digested DNA was separated in a 15 symbol 180 \ f 
“Symbol” \s 12x {20 om 0.7% agarose gel at 1 V/cm overnight in TAE buffer, denatured, 
neutralized, and subsequently transferred to a nylon membrane (Boehringer Mannheim, # 
14172420) by Southern blotting in 10 symbol 180 \ f “Symbol” \ s 12 X } standard saline cit- 
rate (SSC). 

The clones of the three rDNA probes, namely pHr21Aa for 18S, pHrEP for 28S and 
pHr14E3 for INT were originally developed by Dr. M. Muramatsu and were obtained from 
Japanese Cancer Research Bank (JCRB). The positions of each cloned probe in the rDNA ar- 
ray are indicated in Fig. la. For probe labeling and signal detection, we used the Dig DNA La- 
beling and Detection Kit (Boehringer Mannheim, # 1093657), and all the protocols followed 
the manufacture’s specifications. 

2.3 Construction of restriction maps 

Restriction maps of the rDNA spacers for each species were constructed based on hy- 
bridization patterns with three rDNA probes after appropriate single and double digestion. We 
compared our results with the conserved restriction maps of the rDNA coding region constructed 
in previous analyses (Hillis et ai., 1986). Because that the two EcoR I sites on the 3’ ends of 
both 28S and 18S rDNA coding regions (Fig. la) are conserved in vertebrates (Suzuki et al. , 
1990), cleavage site of EcoR I on the 28S downward spacer was determined according to the 
patterns of EcoR [ single digestion after hybridization with 28S probe. Sites of other enzymes 
inside the region were determined subsequently by double digestion of EcoR I with each of the 
given enzymes. For those sites outside this region, double digestions will yield the same pattern 
of a single digestion with EcoR I , and hence were determined by single digestion and by refer- 
ring their conserved cleavage sites in coding region. The same procedure was used to map those 
sites upward the 18S gene. INT probe was used to cleavage sites in the internal spacer region. 
2.4 Phylogenetic analysis . 

The presence and absence of a site on a given position of the restriction map were repre- 
sented by “1” and “0” characters, respectively. The restriction sites among the maps were 
aligned by introducing length variations into the region next to the 3 end of 285 gene (Suzuki 
etal., 1994) to maximize sequence similarity. Thus we obtained a matrix of characters for 
phylogenetic analysis (Table 1). The data set serve to construct a maximum likelihood tree us- 
ing RESTML program in PHYLIP 3.51c package (Felsenstein, 1993) as well as parsimony 
trees using PAUP 3.1.1 software (Swofford, 1991). In all the phylogenetic analyses, the 
lynx was used as an outgroup. 


Table 1 The “0, 1” characters of restriction sites of rDNA spacers in the giant panda and related species 








Animal Characters om 
285 18S INT . 

Black bear G00 1111 10.100 100001 00000000 100001 000d) 10d 1000000 1991 1010100000000004 | 11000101001 \ 
Sun bear ONAT 1111010010000 | 00000000 1000000000 1 1001001000 100110101 000000000 £1 11000101001 1 
Giant panda 0010101 10000 100000 | 0000000000001 1 1001000011001 t 0100101010000 100 1000101000100111 a 
Lesser panda [EE L000001 111011 100000010000 10000000 1001 00000000 OCHONCOO 10 10111001 1000010001 1001 a 
Raccoon 000 10000 1 00000001 100101001 10001 10100006000000000 000001 LOOOL L01 1010000000011 L001 t a 
Lynx G00 1 0000100009001 000101 1010000 100100100001 1010100 COLDODODOO0T 1111 1000000C I TIEIOE Q 





Note: “1” indicates the presence of a site, and “0” the absence of a site. 
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Fig,1 Restriction maps of the DNA coding region (a) and spacer regions (b) of 
the giant panda and related species 
The black boxes are the genes that encode 185, 285 and 5.85 rRNA, respectively. Positions of probes are also 
shown by arrows. With respect to the restriction sites in the spacers, only those nearest the distal end of the gene 
for 185 and 28S rRNA are shown. Asterisks indicate polymorphic sites within the genome of a given animal, 
while the Hindi site downward the horizontal line in the map of the lesser panda is that shows intraspecific varia- 
tion, Abbreviation of species: HB, Asiatic black bear; MB, sun bear; GP, giant panda; LP, lesser panda; RA, 
raccoan; LY, bmx. Abbreviations of restriction sites: A; Stu I, B: BamHI, D: Dm I, E: EoRI, G: 
Bg ll, H: Hind, K: KpnI, L: Bil, P: Pal, S: Sacl, T: BEI, V: Poull, X: Xal. 


3 Results and Discussion 
The rDNA restriction maps of the giant panda and related species are showed in Fig.1. 

Eighteen restriction enzymes, Apa 1, Ava I, BamHI, Bei I, Bll, BEI, Cal, 
Dra I , EcoRI, EcoRV, Hindili, Hpal, KpnI, Pst I ,Pou I ,Sac L ,Stu I ,Xba 1l, 
were used in this study. No maps for Apa I , Awa I , Gal, EcoRV and Hpa I were ob- 
tained due to either one or no cleavage site on the rDNA repeat unit. Therefore the restriction 
maps contain sites for 13 enzymes. By using the INT probe, we detected an additional BstE I 
site in the internal spacer regions of the Asiatic black bear and the sun bear. As regard to the 
intraspecific polymorphism, only one Hind |] site was found differ between the two individuals 
of the lesser pandas in the spacer region near the 18S rRNA gene (Fig.1), The diverged lesser 
panda was not used in phylogenetic analysis. No site was found differ between the two giant 
pandas, hence only one individual is show in Fig.1. The Pou I] pattern of the sun bear always 
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shows two bands when hybridized with 28S probe even after excessive digestion, therefore we 
regard it as a intra-genome heterozygote rather than partial digestion. Both the two sites are 
showed in the map but the variant one was not taken into account for phylogenetic analysis. 

Recent studies showed that the evolutionary rates of both hemoglobin and mitochondrial 
genes among different lineages in the Arctoidea are different (Hashimoto et af., 1993; Zhang 
et al., 1993), therefore when analysing the restriction site data, it is important to employ a 
method which does not necessarily assumes the constancy of evolutionary rate among branches. 
In this study, we used the restriction site maximum likelihood method (RESTML) in 
PHYLIP version 3.51c to construct the phylogenetic tree for the 6 species from the information 
relating to the presence or absence of each restriction site (Table 1). A single tree came out 
with Ln Likelihood = — 173.17 (Fig.2A). In this phylogenetic tree, the Asiatic black bear 
and the sun bear are closely related, then is the giant panda, the lesser panda and the racooon 
form an other branch. 







A: maximum likelihood B: parsimony bootstrap 
Black bear 
100 
Racoon 
97 Sun bear 
Lesser panda 
Giant panda 
Giant panda 
Sun bear Lesser panda 
Black bear Facoon 
Lynx Lynx 


Fig.2 Maximum likelihood tree (A) and parsimony tree (B) for 5 species of Arctoidea 
The lyna is used as en outgroup. The horizontal length of each branch in A is proportional to the estimated di- 
vergence. The figures above some branches in B are the percentage of trees retaining the branch based on 1 000 
bootstrap replicates. 


In addition, we also used exhaustive search option to construct parsimony trees by use of 
PAUP 3.1.1. Three shortest trees were found which have a branch length of 96. These three 
parsimony trees differ in the relation between the lesser panda and the raccoon. In one tree, the 
lesser panda is at the basal branch while in another the raccoon at the basal branch. The third 
parsimony tree which has the same topology with the likelihood tree (Fig.2A) clusters both of 
them in the basal branch. Fig.2B show the consensus tree for the three parsimony tree with the 
support percentage which retains a certain branch in 1 000 bootstrap replicates. This consensus 
‘tree again strongly supports that the giant panda is closer to bear than lesser panda, which is 
consistent with most of the previous research. (O’Brien et al., 1985; Zhang et al., 1993; 
Hashimoto et al., 1993). 

The phylogenetic placement of the lesser panda is another interesting question in Arctoidea 
systematics. O’Brien et af. (1985) suggested that it might be slightly linked to the raccoon. 


342 ae 学 研究 19 4# 





The likelihood tree and one of the three parsimony trees obtained in the present study shows 
that the lesser panda is closer to the racooon than to the bear-giant panda branch (Fig.2). Al- 
though other two parsimony trees failed to cluster the lesser panda and the racooon into one 
clade, their branches are close one another indeed in the trees (not shown). Therefore, we are 
cautiously prone to uhe view that the lesser panda is slightly related to the raccoon. 
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(中 中 国 科 学 院 昆 明 动物 研究 所 细胞 与 分 子 进化 开放 研究 实验 室 昆明 650223) MET: 一 
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摘 要 应 用 rDNA 间隔 区 Southern HERRIRAK ERAT REX 
Ko TAERA, DR, RAR, BRM, BRRR DNA 间隔 区 限制 性 内 切 栈 
图 谱 ， 用 最 大 似 然 法 和 简约 法 构建 它们 的 分 子 系统 树 。 结 果 表 明 大 艰 猫 与 能 具有 较 近 的 闲 
缘 关 系 ， 与 小 熊猫 和 演 熊 的 亲缘 关系 较 远 。 
关键 词 MAN, AMM, DNA 限制 性 图 谱 ， 系 统 发 育 。 六 序 及 | 
中 图 分 类 号 Q959.838 
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“第 十 八 届 国际 遗传 学 大 会 昆明 卫星 会 议 ” 简 介 


“第 十 八 届 国际 中 和 传 学 大 会 昆明 卫星 会 议 ” 于 1998 年 8 月 16 A~18 日 在 云南 省 昆明 市 举行 。 本 次 
会 议 由 中 国 科学 院 昆 明 动 物 研 究 所 与 云南 省 直 传 学 会 联合 主办 ,会 议 主席 由 我 国 著名 植物 学 家 、 中 国 科 
学 院 院 士 吴 征 错 先 生 担 任 ， 秘 书 长 为 中 国 科 学 院 昆 明 动物 研究 所 副 所 长 张 亚 平 研究 员 。 

本 次 会 议 是 一 钦 高 水 平 的 国际 性 学 术 会 议 ， 有 来 自 中 国 、 美 国 、 日 本 、 英 国 、 以 色 列 、 视 大 利 亚 、 
加 拿 大 以 及 中 国 香 港 特 别 行政 区 、 台 湾 省 等 10 余 个 国家 和 地 区 的 知名 专家 和 学 者 约 200 人 参加 。 其 中 
国外 代表 各 余人 ， 国 内 代表 约 170 人 。 会议 共 收 到 论文 捕 要 176 笨 ， 内 容 涉 及 动物 、 植 物 、 微 生物 和 人 
类 遗传 多 样 性 等 诸多 学 科 。 

8 月 16 日 晚 , 会 议 主席 吴征 色 院 士 在 昆明 连 云 宾 迟 主持 了 隆重 的 迎 宾 晚宴， 中 国 科学 院 院士 吴 轴 、 
Refi. RSS, BARS SH Oliver A. Ryder, Ransjit Charkraborty, Rober K. Moizis、Yun-xin Fu, 
吴 仲 仪 和 李 文 雄 等 应 邀 出 席 了 宣 会 。 会 议 于 17 日 在 云南 大 学 开 莫 ， 云 南 省 副 省 长 李 汉 柏 出 席 了 开 基 式 。 
新 华 社 、《 中 国 科 学 报 》、 云 南 电视 台 等 新 闻 单 位 的 记者 也 参加 了 开幕 式 。 开 攻 式 结束 后 ， 图 绕 会 议 主题 
“遗传 学 和 生物 多 样 性 保护 ”， 中 外 学 者 作 了 9 个 大 会 报告 。 著 名 科学 杂志 《Nature》 的 编辑 David Dick- 
son 先生 也 应 邀 作 了 报告 。18 日 ， 与 会 专家 学 者 分 为 人 类 、 动 物 、 植 物 和 微生物 3 个 专题 组 ， 就 种 群 遗 
传 、 进 化 遗传 、 保 护 遗 传 等 内 容 进行 了 广泛 而 深入 地 探讨 。 共 有 名 人 作 了 专题 报告 ， 约 130 人 次 参加 
了 墙报 交流 。 大 会 于 18 日 晚 圆满 结束 。 此 外 ， 围 绕 生 物 多 样 性 资源 保护 ， 于 19 一 22 日 ,会 议 还 组 织 了 
会 后 考 窒 ， 为 中 外 科学 家 充分 了 解 云 南 丰 十 的 生物 资源 、 进 一 步 认 识 云 南 在 开展 分 子 进化 和 遗传 多 样 性 
研究 中 的 独特 地 位 提供 了 更 多 的 交流 机 会 。 

会 议 期 间 和 会 后 ， 数 十 名 中 外 学 者 来 到 昆明 动物 研究 所 参观 和 交流 。 美 国 、 日 本 等 国家 的 几 个 研究 
小 组 与 该 所 科学 家 进 行 了 深入 交谈 ， 商 讨 今后 合作 的 方向 与 有 关 事宜 ， 为 昆明 动物 研究 所 今后 加 强国 际 
间 的 交流 与 合作 拓宽 了 渠道 。 

本 次 会 说 的 举行 ， 对 于 促进 世界 范围 内 的 生物 遗传 多 样 性 保护 与 合作 、 推 动人 类 与 生物 进化 和 遗传 
学 理论 的 发 展 、 以 及 为 我 国保 护 生 物 学 与 进化 生物 学 研究 在 世界 范围 内 赢得 一 席 之 地 ， 都 将 起 到 积极 的 
作用 。 
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